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DISCLAIME R 



Although each program has been tested by its contributor, no 
warranty, express or implied, is made by the contributor or 
1620 USERS Group, as to the accuracy and functioning of the 
program and related program material, nor shall the fact of 
distribution constitute any such warranty, and no responsibility 
is assumed by the contributor or 1620 USERS Group, in con- 
nection therewith. 
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. 1620 USERS GROUP PROGRAM REVIEW AND EVALUATION 

'r.ogram No. Date 

Program Name: 



1. Does the. abstract adequately describe what the program is and 

what it does? Yes No_ 

Comment 

2. Does the program do what the abstract says ? Yes No_ 

Commen t 

3. Is the Description clear, understandable, and adequate? Yes Nc_ 

Comment 

4. Are the Operating Instructions understandable and in sufficient 

detail? Yes No_ 

Comment 

Are the Sense Switch options adequately described (if applicable)? Yes No_ 

Are the mnemonic labels identified or sufficiently 

understandable? Yes No_ 

C omm ent 

5. Does the source program compile satisfactorily (if applicable) ? Yes No_ 

Comment 



6. Does the object program run satisfactorily? Yes No_ 

Comment ; 

7. Number of test cases run 



Are any restrictions as to data, size, range,' etc. covered 

adequately in description? Yes N o 

Comment 

8. Does the Program meet the minimal standards of the 1620 Users 

Group? Yes No_ 

Comment 

9. Please list any suggestions to improve the usefulness of the 
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Commen t _ . 
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THE NUMERICAL CALCULATION OF THE DEFINITE INTEGRAL OF A REAL FUNCTION USING 

SIMPSON'S RULE (Card) 

Donald S. Miller 

Subject Classification 7*0 

Direct Inquiries to : Donald S. Miller 

Applied Mathematics Section 

Management Systems Development Department 

Eastman Kodak Company 

Rochester 4, New York 

Purpose/Description : This program will tabulate using Simpson's Rule, a pre- 
assigned number of definite integrals F(«, c , c^, c ) defined by the 

relation 

x 

F(x, c.j, c 2 , c m ) = J f(x, c^ c 2 , c m )dX 

a 

where l) a, c„ , c„, c are real numbers and 2) f is a real function whose 

evaluation at X, c^, Cg, may be represented in Fortran language by a 

single statement, 

A real number b =» a, and odd subdivision D: a - Xq, X^ + AX^, X^ + 

2AXq, x q + 2nAX^ = b, AK^ « (b-a)/2n, of the interval (a,b) and a range: 

b, b + 2AXq, b + ^AXq, b + 2(k-1 )AX^ for the variable x are associated 

with each F(*,c.,,c , c). 

v 7 r 2 7 ' m 

The numbers m, 2n+1, k, a, b, c.^ c^, • c^ are specified in a 
series of control cards. 

Method : Simpson's Rule is used for approximating all definite integrals. See 
Numerical Mathematical Analysis, Third Edition, by J. B. Scarborough, The Johns 
Hopkins Press, page 152, section h6 and page 17^> section 57 > (c). 
Restrictions /Range : All calculations are subject to the laws of rounding error 
and the precision of the Fortran system- 



if x 



Storage Requirements : N/A 
Equipment Specifications : n/a 

Additional Remarks: This program was written in Fortran language. 
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REGUIAR DATA CARDS 



^-27 



1.23 
I 



+.35124 
II 



0001XXXX 



aooooooooanncoooooooooooooooooooooooooooo 00000000000 00000000000000000000 |||l|||l 

t 2 3 4 5 6 7 I 9 10 1M2 13 14 15 IS 17 10 19 20 21 22 23 24 25 » 27 21 29 30 31 32 33 M 35 39 37 38 39 40 41 42 43 44 45 « 47 40 « ^ 

1 1 1 1 1 | 1 1 1 1 1 11 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 |1 111111111111111111111111 

2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2| 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 1 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2|2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 

3 3 3 3 3 3l 3 3 3 3 3 3 3 3 3 3 3 3 3 3 |3 3 3 3 3 3 3 3 3 3 3 3 3 3 |3 | 3 3 3 3 3 3 3 3 3 3 3 3 ||3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 

4 4-4 4444 444 4444 4444444444444444444444444 W44444444444444|44444444444444444444444 

5 55555555555555555555555555555555555555555555555555 5| 555*55555555555555555555555 
« 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 G 6 G 6 6 8 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 G 6 6 6 G 6 6 6 6 6 6 G 6 6 
7777777777777777777777777777777777777777777777777777777777777777777777777777 SHI 

8 8 8 8 8 8 1 8 8 8 8 8 8 8 8 8 8 8 8 8 8 1 8 8 8 8 8 8 8 8 8 8 8 8 8 8 1 8 8 8 8 8 8 8 8 8 8 8 8 8 8 1 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 

9 9 9 9 9 9 9 9 9 9 9 9 9 9 9 9 9 9 9 9 9 9 9 9 9 9 9 9 9 9 9 9 9 9 9 9 9 9 9 9 9 9 9 9 9 9 9 9 9 9 9 9 9 9 9 9 9 9 9 9 9 9 9 9 9 9 9 9 9 9 9 9 9 9 9 9 9 9 9 9 

1 2 3 4 S 8 7 ( 9 10 11 12 13 14 15 It 17 II IS 20 21 22 23 24 25 26 27 20 21 39 31 32 33 34 35 36 37 3» 39 W 41 42 43 44 45 4S 47 41 49 50 51 52 53 54 55 58 57 51 59 60 G1 S2 83 14 OS 08 17 II 09 70 71 72 73 74 75 76 77 79 70 IS 
UM 5081 



+001.23456; 



oo o||o oo oo oo oo ooooo oooooooooooo oo oo|ooo oiitliooooooooooooooooooo oRJoooo#f|iof||| 

it ii i|i in mi m m ii inn mi nmi|i iiiiii|iiii 11 n 11111111111111111111111111 

2222222|2222222222222222222222222 2 2 2 2 2 |2 2 2 2 2 2 2|2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2f2 2 2 2 
3 3 3 3 3 3|3| 3 3 3 3 3 3 3 3 3 3 3 3|3 333 3 3 3 3 3 3 3 3 3 3§3 3|3 3 3 3 3 3 3 1 3 3 3 1 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 
44 44 44 444|44444 444444444444444444 444 44444 4444 444J44444 4 44444 444 444444444444444 44 

55 55555555|5555555555555555555555555555555555 5555|555 555555555555 555555555555555 
666666666 6 6| 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6| 6666666666 6 6 6 6 6|6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 S 6 6 6 8 

7 7 7 7 7 7 7 7 7 7 7 7| 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7|fff 

8 8 8 8 8 8| 8 8 8 8 8 8| 8 8 8 8 8 8 8| 8 8 8 8 8 8 8 8 8 8 8 8 8 « 1 8 8 8 8 8 8 8 8 8 8 8 8 8 8 |8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 1 8 8 8 8 

Qqqqqqqqqqqqqql 99999 99 99999 9 99999 99 9999 99 999 9 9 999 99 99 9 999 99 999 9999 99 9 9 999 9 999 999 
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Program Description 

Title THE NUMERICAL CALCULATION OF THE DEFINITE INTEGRAL OF A REAL FUNCTION 
USING SIMPSON'S RULE 

Program Input Medium Except for control cards (which will "be described in the 

Program Input section), all input to this program will be punched in regular 

data cards . 

A regular data card is an 80-eolumn, punched card subdivided into four, 
fifteen position fields, columns 1-15, columns 16-30, columns 31-^5 and columns 
J+6-6O, where each sub-field is punched in a right justified manner with a 
decimal number having the form + xxxxx.xxxxxxxx (x represents an arbitrary 
digit). Note that the decimal point is always punched in position seven of the 
corresponding sub-field and that the sign punch may be in any one of the 
positions 1-6. Insignificant zeros need not be punched. The omission of the 
sign punch automatically defines a non-negative number. Any blank sub-field 
will represent zero. Columns 61-80 may be punched in arbitrary fashion and so 
may be used for input sequencing and identification. 

Examples of regular data cards are shown on page 3* 
Program Input For J = 1, 2, • • s let it be required to approximate the definite 
integral 

V X ' °iy C 2J' — ' C mP " fa*' C M> C 2J' "'J c md )dX 

where l) a^, c^, c^y are real numbers, 2) the subdivision D is 

defined by the points = X Q y X QJ + £X QJ , + 2& y X QJ + 2n^X QJ - b^, 

a^ <= b^, ^Xqj - (^j " a ^/ 2n j' 3) the evaluation of the real function f^ at X, 
c^, c 2 y c^ may be expressed in Fortran language by a single statement 
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and k ) the variable x ranges over the k . points b . , b . + 2 AX,- . , b . + 4AX^ , , 
b. +.2(k.-l)AX rt .. 

The input associated with the s functions F 1 , F , •**, F consists of 

the ordered collection C, S , S , S where C represents the Initial Control 

Card and for j = 1, 2, s,S. is the set of all input cards belonging to the 

** J 

3-th function F. . 

(l) The Initial Control Card C is punched according to the following table 



Columns 


Punch 


Description 


1-3 


s 


s is a three 


dig 


Sit 


unsigned number. 


4-6 


m 


i is a three 


dig 


^it 


unsigned number. 



(2) consists of the ordered collection C^ and where C. is the initial 
Problem Control Card and ^ defines the input belonging to the first 



function f ^ . 



(a) The Problem Control Card C^ is punched according to the following 
table 



Columns 


Punch 


Description 


1-3 


2^+1 


2n^+1 is a three digit unsigned number. 


4-6 


k i 


k^ is a three digit unsigned number. 


7-21 




a^ is a thirteen digit, eight decimal 






digit signed number. 


22-36 




b^ is a thirteen digit, eight decimal 






digit signed number. 



(b) The ordered set of coefficients c^, c^, c ^ are punched 

continuously in regular data cards to form the input set ^ . Care 
should be taken to punch the first coefficient c^ in the first 
fifteen position field of the first regular data card in 1 . 
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(3) S 2 comprises the ordered collection and where C Q is the initial 
Problem Control Card and defines the input belonging to the second 
function fg. 

The sets and are defined as in (2), sections (a) and (b). 
(k) The input is completed by punching the sub-inputs S^, S^, S g . 

The ordered collection of sets C, S, , S will then define 

I c. S 

the total input. 

The examples presented at the end of this report should clear up 
any difficulties* 

Evaluation Routine It is necessary for the user of this program to construct 
the so-called Evaluation Routine which is to be inserted in a fixed location 
in the original Fortran source program. 

As before, for j - 1> 2 P s let it be required to approximate 

the definite integral 

x 

F i {x ' °U' °2J' °mj ) * / f J (X ' C n' °2J' "•' °mj )aX • 

The Evaluation Routine is written in Fortran language and will compute 
the value of the function f^ at X, c^, c^y • c^ • In writing this program 
one assumes the argument X to be stored at "Fortran" location X while the 
coefficients c« c ., • c . are stored in the so-called C-array. 

2j Sl_ 

For simplicity suppose that the coefficients c^y c^y c^ are 

stored at locations C(l), C(2), C(M) respectively. 

The user must evaluate the given function f . using the coefficients 
stored in the C-array and the argument stored at X. 
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The resultant functional value must be stored at "Fortran" location Y. 

Tne ordered set of Fortran cards defined by the Evaluation Routine is 
then filed in the original Fortran source program immediately following the card 
with identification (See, Fortran Program Listing at the end of this 

report . ) 

The augmented Fortran source deck may then be used to compile the 
desired object program. After compilation^ the user should remove the Evaluation 
Routine from the source program, which will then be ready to accept a different 
Evaluation Routine. 

Restrictions Refer to the notation of Program Input section. 

The number s must satisfy the inequality 1 25 s 2 999 • 
The number of coefficients m belonging to the integrand f . } 1 2f j 2f s 
must presently satisfy the inequality 1 — m —50. The size of m is controlled 
by the single dimension statement with identification 01020C50. Consequently 
the range of m may be decreased or increased depending on the amount of computer 
memory available. 

For 3 - 1, 2, s, the number 2n^+1 of points of subdivision of the 

interval (a., b.) and the number k. of entries in the tabulation of the definite 

y y j 

integral F.(», c 1 , c c ., c .) must satisfy the inequality 1 — 2n.+1, k. ~ 

999- 

All calculations are subject to the laws of rounding error. 



The identification is punched in columns 73-80 of each Fortran source card. 
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Program Output 

Switches All program switches may "be set arbitrarily. 
Content The notation of Program Input section will "be used. 
For 3 * 1, 2, s, this program will record 

(a) the number 2n.+1 of points of subdivision into which the interval (a., b.) 
is subdivided, 

(b) the number 3^ of entries in the tabulation of the definite integral function 

?(', C^, C 2 y C^), 

(c) the interval limits a., b., 

y y 

(d) the function coefficients c,., c,-,., c 

v ' iy zy 9 my 

(e) the tabulation of the definite integral function F(*, c.^, °2y c mj^* 

Note that the columns headed by X, AREA and E respectively define the 
value of the independent variable x, the corresponding "Simpson's Rule" approxi- 
mation to the definite integral £^( x > °iy C 2$' C m^ an< * ^ e error generated 
in using such an approximation. 

Error Messages The notation of Program Input section will be used. 

The message THERE IS AN INVALID PARAMETER will be recorded if for any 
3 such that 1 — j — s, one of the following error conditions occurs: 

(a) the total number of points of subdivision N in subdivision D is even, 

(b) the total number of points of subdivision N. in subdivision D is less 

j 

than 3> 

(c) the total number of entries in the tabulation of the j-th definite 
integral function c^, c^y c^) is negative. 

Whenever one of the above error conditions arises, the given program 
will either proceed with the next numerical integration or will stop at end-of- 
Job. 
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Numeric Example 

1 . Evaluate numerically the definite integral 

x 



J x 2 Vi+6x dX 



for x = 1, 1.04, 1.08, 1.76 . 

2. Tabulate the function defined by the integral 

x 



over the range 5 — x — 6 . 
3. Compute numerically the definite integral 



k. Compute numerically the definite integral 



To apply the given program to examples 1 and 2 define the function f 
by the relation 




27T 




7T 




2 



F(X, Cl , c 2 ) X 2 V 



C+c^X 



A direct verification shows that 



x 



x 












x 



X 
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Let it be assumed that the coefficients c^ and c^ are stored in the 
C-array at locations C(l) and C(2) respectively. 

Evaluation Routine 1 for evaluation f at X, c^ and is shown on 

page 18. 

To apply the given program to examples 3 and k define the function g 
by the relation 

A direct verification shews that 

?1X 2jX 
| r- 32 = |g(xAl,1,0)dX , 

6 X + sin 9 

n ix 

2 cos 9 



Let it be assumed that the coefficients c 1 , e , c x and Ci, are stored 
in the C-array at locations C(1), C(2), C(3) and 0(4) respectively. 

Evaluation Routine 2 for evaluating the function g at X, c^, c^, 
and c^ is shown on page 19 ♦ 

Evaluation Routine 1 must then be inserted in the Fortran source 
program immediately after the card with identification (columns 73-60) 3011 C50, 
After filing this evaluation routine, the resultant source program may then be 
compiled using the Fortran processor to give Object Program 1 . A similar 
operation performed on Evaluation Routine 2 will generate Object Program 2 , 

Input 1 for examples 1 and 2 will be used with Object Program 1 , 
s and m will both have value 2. 

^^^^^ 
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For example 1, set « 0, b 1 ■ 1, 2^+1 ■ 51 , * 20 and B^^ ■ 0>l) 
For example 2, let 8g » 1, bg « 5> 2ng+1 « 151, kg = 25 and Sg * (7,9) 
Input 2 for examples 3 and k will be used with Object Program 2. 
s and m will have values 2 and k respectively. 

For example 3, set a 1 « 0, b 1 = 27T, 2^+1 = 75, ^ = 1 and Sj = 

For example k, set a « 0, b = 7T, 2n +1 = 125, = 1 and S . = 

d. d. d d d\ 

\2> > * 2 

The tabulations for the integrals 

x x 

+6X dX and J X 2 V7+9X dX 
1 

are shown in Output 1, Problems 1 and 2 respectively. 
The values of the integrals 



27T 7T 

J ^ + sin 9 i 



J | + sin 9 ^ | - cos e 

are shown in Output 2, Problems 1 and 2 respectively. 



b„ - a. 



* 11 

Note that Ax, - .02 and since Ax. ■ -~ , then 2n„+1 » 51. 

1 1 2n^ 7 1 



c 



o 



INPUT 1 



c 



o 
o o 
o o 



X X 
X X 



cvi m 
o o 



o 
o 
o 
c 
o 
c 
o 
o 



c 

Oi 
O! 

c 

<c 
c 
o 
o 



— ; If i 



I 

• \ • 

o -«jo o 

o ,o 
c 



o 
o 
o 
o 



'Oi 

;o: 
oi 

•I I •! 



cvj o «-< in r- 

o| CVJi I CVJ 
O! 0| i O 
(VI -«j 
O; If) 

o! o 



in 



o 



c 



OOTPOT 1 




I 



C 



O - 



NUMBER OF COEFFI C TENTS EQUALS 



"PROBLEM NUMBER 



N K A B 



51 20 .00000000 1.00000000 



function coefficients 



COEFFICIENT NUMBER COEFFICIENT 



1 1.00000000 



1 .00000000 





INTEGRAL 


FUNCTION 






1 

* 




X 


AREA 


E 


1_ 




1.00000000 


,44 C 24 1 82 


•00000000 






1 • 04000COO 


.49940039 


•00000000 






1 • 08000000 


.56391 849 


.00000000 






* 1.12000000 


.63406824 


•ccooooco 






1.16000000 


•71 012593 


•00000000 






1 .20000000 


.792371 93 


•00000000 






1 .24000000 


.88 1 09052 


•ooooooco 






1 .28000000 


• 976570 1 2 


-.00000000 






1 .32000000 


1 • 0791 0290 


-.00000000 






1 .36000000 


1 . 1 8898490 


-.00000000 






1 .40000000 


1.30651620 


.00000000 






1.44000000 


1 .43200030 


-.00000000 






1 .48000CC0 


1 .56574460 


-•oooooooc 


— g o 




1.52000000 


1 .70806020 


-.00000000 


— £o 




1 • 56000000 


1.85926180 


-•oooooooo 






1 • 60000000 


2.0196676 


•00000000 







1.64000000 


2 • 1 8959970 


-.oooooooo 


OC 




1 .68000000 


2.369383 30 " ■ 


-.oooooooo 




1 .72000000 


2.55934720 


-.oooooooo 


F= — h 




1 .76000000 


2 .759324 00 


-.oooooooo 


Orn 

- -< - 











PROBLEM NUMBER 2 

- ■ - - ■- -• ' \ 

N ~~ K ~ A B " , 

151 25 1 .00000060 5.00000000 " 



Function' Coefficients 

coeff icient"' number' cofff i c i ent 

__ _ 7.00000000 

9.00 0000 



INTEGRAL FUNCTION 



X 


AREA 


E 




5. 00000000 


262.98258000 


-.00000008 




5.05333330 


272.72274000 


-.0000001 1 




5. 1 0666660 


282.71626000 


-.00000002 




5.15999990 


292.96706000 


•00000002 




5.21333320. 


303.47907000 


-.00000002 




5.26666650 


314.25626000 


•00000005 


1 


5.31999980 


325.30256000 


.00000008 


-f=- 


5.37333310 


336.62 1970 00 


•00000005 


1 


5.42666640 


348.21849000 


•00000002 




5.47999970 


360.09614000 


• 00000002 




5.53333300 


372.25898000 


-•00000005 




5.58666630 


384 .71 1 05000 


-.00000002 




5.63999960 


397.45644000 


-•00000002 




5.69333290 


4 1 0.49927000 


•00000002 




5. 74666620 


423.84363000 


-.00000017 




5.7999995C 


437.49368000 


-.00000017 




5.85333280 


451 .453550 


-.00000002 




5.90666610 


465.72744000 


• 00000005 




5.95999940 


480.31954000 


.00000026 




6.01333270 


495 .23405000 


• 000000 1 1 




6.06666600 


510.4751^000 


•00000014 




6. 1 1999930 


526.04726000 ' 


. 00000020 




6.17333260 


54 1 .95448000 


.00000038 




6.22666590 


558.20 108000 


-.00000029 




6.27999920 


~ 374.79150 000 " 


" .00000071 
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OUTPUT 2 











NUMBER OF COEFFTCTElTrs~XQUAT_S 4 



PROBLtM NUMBLR 



N K A B 



75 1 .00000000 6.28318530 



FUNCTION COEFFICIENTS 
COEFFICIENT NUMBER COEFFIC IENT 



1 1.25000000 



1 .00000000 



1 .00000000 



00000000 



I 



INTEGRAL FUNCTION 




ON 
■ 




X AREA 


E 





6.28318530 8.37757770 -.00000002 



PROBLEM NUMBER 



^ K A B 

125 " I .00000000 3. 141 59260 



FUNCTION COEFFICIENTS 
COEFFICIENT NUMBER COEFFICIENT 



1 1 .50000000 

~2 -1 .00000000 



"3 -1 .00000000 

"5 1 .57079630 



INTEGRAL HUNLT10N 



C) 

o 
c 
c 
c 

o 
o 
o 



iO 

cc 

(VI 

o 

03 



o 

v0 

|CVJ 

a- 
in 



EVALUATION ROUTINE LISTING 1 



STAT FO RTRAN 

NO • STATEMENT 



C 



EVALUATION ROUTINE 1 



SER PRO 
NO* NO. 



400 I OC 50 
40020C50 



00 
-I — 



EVALUATION ROUTINE LISTING 2 



G 



© M 



_STAT_ FORTRAN SER PRO 

NO. STATEMENT " " NO. NO. 



EVALUATION ROUTINE 2 40010C50 



Y=l ./(C<~1 >+C(2)*SlNF<C(3)*X+C<4) ) ) ~ 40Q2QC5Q 



i 



VO 

-»- 



m -i 



o 



o 



FORTRAN PROGRAM LISTING 



.J" 



W 


STAT 
NO. 


... _ _ __ . FORTRAN 

STATEMENT * 




SER PRO 0^ 

NO. NO. . V^._ 

- - 


- - 

— - 


_ _ 




- — - — 




- - 


- - 


c 


-~ - 


- _ ... - ... ... - — - - — . — - — - - . 

AREA UNDER CURVE USING SIMPSONS RULE 




000.1 0C50 


• 


c 




DIMENSION STATEMENTS 




01010C50 








DIMENSION C(50)»S(6) 




01 020C50 





c 




FORMAT STATEMENTS 





020 1 0C50 


— : 




1 


F0RMATC2I3) 




02020C50 


■ 




2 


FORMAT (1H ///38H NUMBER OF COEFFICIENTS EQUALS 14) 




02030C50 







3 


FORMAT (1H ///30H PROBLEM NUMBER 14) 




02040C50 






4 


FORMAT ( 2 1 3 • 2F 1 5 • 8 ) 




02050C50 






5 


FORMAT (1H /49H N K A 


B) 


02060C50 






6 


FORMAT ( 5H I4»4H I4t5H F15.8»5H F15.8) 




02070C50 






7 


F0RMAT(4Fl5.8) 




02080C50 






8 


FORMAT ( 1 H /35H FUNCTION COEFFICIENTS) 




02090C50 






9 


FORMAT (44H COEFFICIENT NUMBER COEFFICIENT) 




02 1 00C50 






1 


FORMAT (13H I3»15H F15.8) 




02 1 1 0C50 






1 1 


FORMAT (1H ///33H INTEGRAL FUNCTION) 




02120C50 






12 


FORMAT (1H /49H X AREA 


E ) 


02 1 30C50 






13 


FORMAT (F19.8»F19.8»F19«8) 




02 1 40C50 






1 4 


FORMAT* 1H /39H THERE IS AN INVALID PARAMETER) 




02 1 50C50 


1 

to 


c 




READ IN PROCESS NUMBER CONTROL CARD 




040 1 0C50 




1 — 




999 


READ 1«LN0P»LCN0 




04020C50 








PUNCH COEFr ICItNT NUMBER 




050 1 0C50 








PUNCH 2»LCN0 




05020C50 




L. 




COUNT PROBLtWS 




060 1 0C50 








DO cJZTj LPNO= 1 t LNOP 




06020C50 




L. 




HUNCH PKUBLtn NU^btK 


- 


070 1 0C50 








PUNCH 3 % LPNU 


- 


07020C50 








KtAU IN PROuLtiM CONTROL CARD 




080 1 0C50 








READ 4»NiK»AiB 




08020C50 




C 




PUNCH PARAMcTtRS 




090 1 0C50 








PUNCH 5 




09020C50 








PUNtn o»N»ls»A«D 




09030C50 




L 




READ IN COEFFICIENTS 




1001 0C50 








DO £-0 1- I »LtNO<4 




1 0020C50 






*5 f\ r\ 

200 


KcAU T*v-( I ) 1 + 1 ) »t( »L( 1+3) 




10030C50 








Di i Nlf~ ui rnccc T f~ T Fmtc 
r-'UiNv-rl CUur r l u l CN I o 


- - 


1101 0C50 









PUNCH 8 


— 


1 1 020C50 


.. 






rUNtn V 





1 1 030C50 


. - 






UU cOl I-l«Lv-NU 




1 1 040C50 







20 1 


PUNCH 10»I»C(I) 





1 1 050C50 




c 




ERROR ANALYSIS 













ODD CHECK N 




13010C50 





STAT 
NO. 


FORTRAN 
STATEMENT 




SER PfkC 
NO* NO. 









IF (N-( N/2 ) ) 223 « 223 * 202 




1 3020C50 





C 


_ _ 


N IS ODD .NON-NEGATIVE CHECK N-3 




14010C50 






202 


IF(N-3)223«226,203 




i4p2qcj:o 




C 




N EXCEEDS 3. SET N=3 BIPASS SWITCH AND I NC=-3 




15010C50 


; 




203 


JTH1=1 




15020C50 








lNC=-3 




15030C50 




C 




NON-NEGATIVE CHECK N-5 


_ 


1601 0C50 









IF (N-5) 226*224 • 204 




16020C50 




C 




N EXCEEDS 5. SET N=5 BIPASS SWITCH 




17010C50 






204 


JTH2=*1 




1.7020C50 




C 




POSITIVE CHECK K 




1801 0C50 






225 


IF(K)223.223,205 




18020C50 




C 




K IS POSITIVE. CALCULATE H 




19010C50 






205 


AN M 1 = N- 1 




19020C50 








H=(B-A)/AN M 1 




19030C50 




C 




POSITIVE TEST H 




20010C50 








IF(H)223.223,206 




20020C50 




C 




H IS POSI TIVE. PUNCH TABLE HEADINGS 




21 010C50 


rb- 




206 


PUNCH l l 




21 020C50 


1 — 






PUNCH 12 




21 030C50 




C 




INITIALIZATION 




22010C50 








JB= 1 




22020C50 








JE=N+2 




22030C50 








X=A-H 




22040C50 








JTH=1 




22050C50 








YOD=0.0 




22060C50 








YEV=0.0 




22070C50 








JTH3=1 




22080C50 




c 




COUNT NUMBER OF AREA VALUES TO BE COMPUTED 




23010C50 








DO 211 I = 1 ,K 




23020C50 




c 




COUNT NUMBER OF DEPENDENT FUNCTIONAL VALUES TO BE 


COMPUTED 


24010C50 








DO 208 J=JB,JE 




24020C50 




c 




COMPUTE L 




25010C50 








L= J-JE+6 




25020C50 




c 




COMPUTE FUNCTIONAL VALUE AT X AND STORE AT Y • NOTE 


THAT THE 


30010C50 





c 





FUNCTION COEFFICIENTS ARE STORED IN C-ARRAY 




3001 1C50 




_ _ _ 




BIPASS SWITCH 

GO T0C207, 218*222) ,JTH 




6001 0C50 
60020C50 




c 

6 


207 


STORE Y( 1 > .Y<2 > AND YC3) AT S < 1 )'» S ( 2 > AND S ( 3 ) 
S( J)=Y 
3 TEST J 

■ " ■' " - © 


- — - - 


61 010C50 
61 020C50 
62010C50 

O 






Ci 

CT AT CADTD A M " - • 



_SXAT _ FORTRAN ^ SER PRO 

NO. STATEMENT NO. NO. 







I F ( J-3 ) 208 *2 1 6 * 223 


62020C50 


: 


c 


- - 


J IS. LESS THAN 3. INCREMENT X BY H 


630 1 0C50 






208 


X = X+H 


63020C50 


- - 


c 




COMPUTE ERROR AND AREA 


640 1 0C50 




c 




N=3 BIPASS SWITCH 


650 1 0C50 









GO TO (209 « 2 1 4 ) » JTH 1 


65020C50 





f 




N~5 BIPASS SWITCH 


ddO 1 0Lr>U 









AS V AN ^ A*% « AN AN 4 AN % a * J 

GO TO ( 2 1 ♦ 2 1 2 ) • JTH2 


66020C50 




c 




CALCULATE ERROR FOR N GREATER THAN 5 


670 1 0C50 







210 


ERR0R=S( 1 )+S(6)~4.*(S(2)+S(5) ) +7.* ( S ( 3 ) +S ( 4 ) ) -8 • * ( Y0D-YEV ) 


67020C50 








ERROR= ( — H/90 • ) *ERR0R 


67030C50 





c 




CALCULATE AREA FOR N GREATER THAN 3 


680 1 0C50 


- 




213 


AREA= (H/3. )*(S(2)+S(5)+4.*(S(3) +YEV+S ( 4 ) ) +2 .*Y0D ) 


68020C50 




c 




PUNCH . I NDEPENDENT VAR I ABLE % AREA AND ERROR 


69010C50 






215 


PUNCH 13»B»AREA , ERROR 


69020C50 




c 




INCREMENT B BY 2H 


700 1 0C50 








B=B+2.*H 


70020C50 




c 




UPDATE Y0D AND YEV 


7101 0C50 


1 






YOO=:YOD+S ( 5 ) 


71 020C50 


ro 

1 — 






YEV=YEV+S(4) 


71 030C50 




c 




1— N fw t**\ 1 A A^> r"~ v/ i kl A T* *•*#./«% r~> v/ v a A fc. 1 • « k A "f A"» * A* » 

REPLACE Y(N— 1 ) AT S(4) BY Y(N+1 ) AT S(6) 


720 1 0C50 








5(4) =S ( 6 ) 


72020C50 




v» 




1 inn A T C 1 Q A Kin 1 FT 

UrUATh AND Jt 


730 1 0C50 








ID — ICj. 1 


73020C50 






O 1 1 

2 11 




73030C50 


— _ 

— 






GO TO 2.2. r 


73040C50 




c 




bt 1 IN = tD DlrAbo owl 1 Ln»LALLULATt EKKUK r OK IN = 5 


740 1 0C50 






212 


JTnt 2 1 


74020C50 








tKKuK= ( -n/yt) . J*(b( 1 ;+o(o)— 4.*(b(2)+b(5; )+f.*(S(3)+S(4) ) — 8 »*YUD ) 


74030C50 


: — 






v^O I O c 1 J 


74040C50 











bb 1 !N = 3 dIKAbb owl I Ln»LALLULATt tKKUK Ainu AKtA rOK N = J 


7o0 1 0^-50 







214 


IT" Ul 1 — 1 

J I HI - 1 


■7C f\ *5 f*\ cz r\ 








tKKUK = \ ~rlx yu • )*ibl 1 )+blo)*'*«*lb((;)Tb(b) )+o.*b(J / ) 


75O30C50 




- - — - 




— — 






— 


c 




SET S ( 4 ) =0 


0C50 








S(4)=0.0 


76020C50 








GO TO 215 


76030C50 




c 




J=3.SET BIPASS SWITCH. TEST L+INC 


7701 0C50 






216 


JTH = 2 


77020C50 






220 


IF (L+INC ) 208,217,223 


77030C50 




c 




L+INC=0.SET BIPASS SWITCH 


' 78010C50 





FORTRAN 
STATEMENT 



SER PRO 
NO. NO. 





217 


JTH = 3 


78020C50 ' 








GO TO 208 


78030 C50 




c 




YOD-YEV B I PASS SWITCH 


7901 0C50 _____ 






218 


GO T0<219«221 ) » JTH3 


79020C50 




c 




SET YOD-YEV BIPASS SW I TCH . UPDATE YOD 


80010C50 






219 


JTH3=2 


80020C50 








Y0D=Y0D+Y 


80030C50 








GO TO 220 


80040C50 




c 




SET YOD-YEV BIPASS SWI TCH. UPDATE YEV 


81010C50 






221 


JTH3=1 


81020C50 








YEV=YEV+Y 


81030C50 








GO TO 208 


81040C50 




c 




STORE Y(N)«Y(N+l) AND Y(N+2) AT S(4)»S(5) AND S(6) 


82010C50 






222 


S(L)=Y 


82020C50 








GO TO 208 


82030C50 




c 




H IS NOT POSITIVE. THERE IS AN INVALID PARAMETER 


83010C50 






223 


PUNCH 14 


83020C50 








GO TO 227 


83030C50 




c 




N*5.SET N=5 BIPASS SWITCH 


84010C50 


| 




224 


JTH2=2 


84020C50 








GO TO 225 


84030C50 




c 




N=3.SET N=3 BIPASS SWITCH AND I NC=-4 


85010C50 






226 


JTH1=2 


85020C50 








INC=-4 


85030C50 








GO TO 224 


85040C50 




c 




CONTINUE STATEMENT 


86010C50 






227 


CONTINUE 


86020C50 




c 




END OF JOB 


87010C50 








PAUSE 


87020C50 








GO TO 999 


87030C50 








END 


87040C50 





-24- 



General Remarks 

1 . Definite integrals may be approximated without using any coefficients in 
the C-array. For example, let it be required to tabulate 




for x = 1, 1.02, 1.0*t, 1.06, 1.08, 1.10. 

The Evaluation Routine can be defined by the single Fortran 
statement 

Y = EXPF(-X*X) . 

Referring to the notation of the Program Input section, the input for 
the corresponding object program would consist of (a) a Process Number 
Control Card C with m = s = 001 , (b) a Problem Control Card with 
2x^+1 = 51, k = 6, a 1 = 0, b <| = 1 and (c) a single blank card. 
2. (a) The following referents must not be used in programming the evaluation 
routine 

1-1 h inclusive, 
200-227 inclusive, 
999 . 

(b) The following variables must not be used in programming the evaluation 
routine 

A, ANMI, AREA, B, C, ERROR, H, I, INC, J, JB, JE, JTH, JTH1 , 
JTH2, JTH3, K, L, LCN0, LNOP, LFNO, N, S, YEV, YOD. 
5. If after tabulating a series of definite integrals, it is found necessary 
to make some changes or to investigate a new set of definite integrals (of 
the prescribed type), it is not required to re-read the object program into 
memory. Merely feed the new input into the card reader, press the console 
start key and then the card reader, reader-start key. 
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